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ABSTRACT 

We study the orbits in the MOND theory within a dwarf galaxy of mass 
M,j ~ 10 8 Mq at a distance of ~ lOOkpc from a neighboring galaxy of mass 



5 x 10 -Mq, such as ours. It is assumed that a second mass m « is 
gravitationally bound to by a previously calculated potential for the MOND 
theory. This potential is obtained for a free falling mass in a constant external 
gravitational acceleration field V0 9 . The numerical technique of surfaces of section is 
used to study the stability of the phase-space orbits in the dwarf galaxy. Equatorial 
orbits with sufficiently small eccentricities e < 0.65 are found to be stable with respect 
to small changes in the initial conditions. (The equatorial plane is perpendicular to 
the direction of V</> 9 , which is along the line joining Md and M g .) For decreasing 
values of the conserved component of the angular momentum, in the direction of V(/> 9 , 
equatorial stability is lost. 

Subject headings: cosmology — gravitation — galaxies: dwarf. 



1. Introduction 

The inconsistency between luminous and dynamical mass measurements is well known. An 
alternative way to explain this inconsistency without dark matter, is by Modified Nonrelativistic 



Dynamics (MOND) for galaxy and galactic systems (Milgrom 1983a, 1983b , 1983c) . This theory 



was later developed into an alternative theory of gravitation by Bekenstein fc Milgrom (1984), 



which we hereafter refer to as the Bekenstein-Milgrom (BM) theory. MOND agrees with the Tully 



fc Fisher (1977) and Faber fc Jackson (1976) laws, and explains the dynamics in elliptical and 
spiral galaxies, without the need for dark matter. It satisfies conservation of energy, momentum, 



angular momentum, as well as the weak equivalence principle (Bekenstein Sz Milgrom 1984). 
Recent studies of rotation disks of Low Surface Brightness (LSB) Galaxies by [Salcedo "& 



Gamez (1999)], |de Block fc McGaugh (1998)| and |McGaugh fc de Block (1998)| are in agreement 



with MOND. 
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Departure of MOND from Newtonian theory is not connected with a distance scale. In 
MOND, deviations from Newtonian theory become significant for very small accelerations. The 
main argument presented is that Newtonian gravitation has not been tested for very weak fields 
and it could be that the theory is not valid in this regime. Other theories modifying Newtonian 
theory were reviewed by BekensteTn] (1987). 



Stability of galactic disks in BM theory was first studied in a WKB approximation by Milgrom' 



(1989) . The WKB scheme deals with perturbation wavelengths A = 2ir/k, much smaller than the 



involved distances p on the disk (\kp\ >> 1). In Newtonian gravitation, the presence of a dark 



halo is important in stabilizing the disk against violent bar formation (Ostriker fc Peebles 1973) 
Recent N-body simulations for the BM theory, indicate that disks are more stable in MOND than 
in Newtonian dynamics with dark halos ( Brada fc Milgrom| 1999). 



We obtain the conservative Hamiltonian that describes the motion of a particle in a previous 
calculated potential for the BM theory in Section |2[ This Section also includes a summary of the 
stability theory for phase-space orbits. Orbits in a dwarf galaxy are discussed in Section |3[ Our 
conclusions are presented in Section ||. 



2. Theory 

In the BM theory, the Poisson equation for determining the gravitational potential is modified 

to 

V. (p (^^) Vy>) = AttGp, (1) 

where ao = 2 x 10 _8 cm/s 2 is a constant with the dimension of acceleration, set so as to agree with 
the Tully-Fisher law (Milgrom 1983b|) , the function p(x) (where x = |V</?|/ao) obeys < p(x) < 1, 



with lim^^o ^( x ) = x an d G is the Newtonian gravitational constant. With this particular value 
of ao, non-Newtonian effects due to the Solar gravitational field are expected to begin only beyond 
the Oort cloud |(Oort 1963). 



Let us assume the existence of a constant external gravitational acceleration V(/> 9 due to some 
source (e.g., a neighboring massive galaxy, such as ours) and a free falling sphere of mass M d (e.g., 
the center of a dwarf galaxy). Sufficiently close to M d so that p(x) 1, the solution of equation 
(|I]) yields the well known Newtonian potential 

GM d 

ip « . 

r 

When the solution tp of equation (|l]) implies gravitational accelerations much bigger than 
the external gravitational acceleration (\Vp\ » |V0 ff |) but \V<p\ < ao, spherical symmetry can 
be assumed and the potential is obtained from equation (||) in the limit p(x) » x, using Gauss's 
theorem 

ip^ VGM d a logr. (2) 
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In the opposite limit, when |V</>| << |V</> 9 |, Bekenstein fc Milgrom| (1984) have calculated the 



approximate potential of equation (|T|) in cylindrical coordinates 



^(1_ Q )+p 2 ' 



where z is the direction of V<^ s , /i 9 = /i(| V(/> 5 |/ao), £ 9 = dln(^)/dln(x) U=|V0 9 |/a o an d 

a 9 = L g /(1 + In the asymptotic Newtonian limit |V</> 9 | >> ao, /x 9 — > 1, L g — > 0, a 9 — ► 0. In 

the non-Newtonian limit, which we call the MOND limit, |V</> 9 | << ao, 

fig(x) ~> X (4) 

L 9 (x) - 1 (5) 
a 9 -» 1/2. (6) 

It is to be noted that equation (||) is not a consequence of tidal effects, since the external field is 
constant ( Bekenstein Milgrom| 1984). 



We study the orbit of a very small mass m « such that it does not modify the potential 
in equation (||) and consider the motion of m, bounded by this potential. This is a central force 
type problem, so that conservation of linear momentum implies that the coordinate dependence 
of the Hamiltonian is in the relative distance between and m. According to equation (||), and 
considering the azimuthal symmetry of the potential, the dynamics of the — m system are 
governed by the conservative Hamiltonian 



" Pl+Pl + Vj- 1 A (7) 



2// 



A = (p g y/l + Lg)- 1 GM li m, (8) 

where [x r = M^m/ {M^ + m) ~ m is the reduced mass and l z is the angular momentum component 
in the direction of V0 g (the direction joining the neighboring massive galaxy and the dwarf 
galaxy). Thus, the three dimensional motion in a axisymmetric potential is reduced to the 
motion in a plane. In equations ([?]) and (11), cartesian coordinates (p,z) are used to describe this 
(non-uniformly) rotating plane, which is often called the meridional plane. 

In the MOND limit, when |V0 S | << ao, equation (@) with the condition |V^| << |V0 g |, 
yields a minimum distance between and m: 



r -J GMd a ° (9) 

where the values for L g and p g are obtained from equations (Q)- (||). We note that the forces in 
the p and z directions are not equal. 
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The Hamiltonian (|7|) has an elliptic equilibrium position at z = 0, p = /^/(/i r A), which is 
called the guiding center. Expanding the Hamiltonian in a power series about the guiding center, 
we have 



H = H (I)+H 1 (I, 
H (I) = u z I z +iv I, 



p^p 



(10) 
(11) 



This expansion is called the epicycle approximation and is appropriate in the neighborhood of the 
equilibrium. Hq is the unperturbed motion and describes the epicycles around the guiding center; 



the two frequencies uj p and uj z given in equation (|ll|) , are called the epicycle frequency and the 
vertical frequency, respectively. In phase space, the conservation integrals define a 2-dimensional 
torus. When j z uj z + j p uj p = 0, j = \ j z \ + \j p \ for integers j Z)P , the linearized frequencies are said to 
satisfy a resonance of order j. The corresponding torus is called a resonant torus. Prom equations 



(|6|) and (|TT]) , the allowed linearized resonances are restricted to a g < 1/2. 

For a resonance of order j, it is possible to put the Hamiltonian into a Birkhoff normal form 
of degree j by a canonical transformation 

H = K (J) + K l (J,<!>), 



where Kq is a polynomial of degree j in the new actions J and $ are the new ang les |( Arnold 1989)1 
Kq is a good approximation when the perturbation K\ is sufficiently small. The new frequencies 
are, then, 

_ dK (J) 



(12) 



and are functions of the new amplitudes J% , a phenomenon first discovered by [Lindstedt (1882)| in 
connection with nonlinear oscillators. 



Nondegenerescence is defined as the nonvanishing of the Hessian determinant 

dvi 



det 



dJ k 



(13) 



There exists a theorem of Kolmogorov (1954) on the behavior of a nonresonant torus under a 
small perturbation K\ of a nondegenerate Hamiltonian Kq (Arnold 1989), which was subsequently 
proved by [Arnold (1961) for Hamiltonian systems and by [Moser (1962) for area-preserving maps 



(see also Arnold 1963 and Moser 1971 ). The theorem, known as KAM in recognition of their work, 
states the existence of invariant tori, densely filled with phase space curves winding around them, 
which are conditionally-periodic with the number of independent frequencies equal to the number 
of degrees of freedom. This theorem is valid if the following condition holds 



\{y.j)\>C\j\\ 



(14) 
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In equation ( |i~4| ) C is dependent on the magnitude of the perturbation and r, on the number of 
degrees of freedom. The existence of an invariant, or KAM, torus surrounding a periodic orbit, 
characterizes the stability. Orbits whose frequencies satisfy equation (|14|), occupy most of the 



volume in phase space for sufficiently small C (Moser 1971) 



The conserved l z is an isolating integral for the full three degrees of freedom system; the two 
degrees of freedom system described by equation (|7|), is a subspace of constant l z of the original 



three degrees of freedom system. A surface of section can be constructed in the spirit of Henon 



fc Heiles| (1964). In a two de grees of freedom system, for instance equation (0), the dimension 



of the phase space is four. Since the energy is an isolating integral, motion is constrained to a 
three dimensional constant energy surface, for example, p, z and p z . Successive intersections of 
the trajectory with the plane p = const, and p p > 0, described by the set of points z and p z (or 
z = dz/dt), is called a surface of section. Each intersection of the orbit is a point in this plane 
and the passage of one point to the next can be considered as a mapping. After an infinite time 
interval, the points corresponding to a unique orbit fill up a whole area of the surface of section, 
in the absence of further isolating integrals. 

Motion is constrained to the intersection of the surfaces of constant H and I in the presence 
of a third isolating integral. Consider an orbit whose initial conditions lie in the phase space region 
of conserved /. Its points in the surface of section form a smooth curve. 

According to the Poincare-Birkhoff theorem for resonant tori, after the inclusion of the 
perturbation K\, we have an even number of periodic orbits in the vicinity of a stable periodic 
orbit. Periodic orbits are exact resonances of the nonlinear problem. The corresponding surface 
of section is a fixed point, surrounded by an even number of other fixed points. Half of them are 
stable, while the other half are unstable. Stable points are surrounded by closed invariant curves. 
Unstable points are connected by separatrices, forming a structure called a chain of islands. Chaos 



is present in the vicinity of the separatrices. (For further details, see Lichtenberg Sz Lieberman 
(1991)| and |Dunby (197i].) 



3. Orbits in a Dwarf Galaxy 

In the following, we describe the motion of a particle m (e.g. a globular cluster) around a 
spheroid of mass (the nucleus of a dwarf galaxy) with M& » m, as shown in Figure [l|. We 
take 

M d = 1O 8 M , 
m = 1O 5 M . (15) 

We assume that the Mj — m system is gravitationally bounded to a nearby massive galaxy 
M g = 5 x 10 11 Mq (such as the Milky Way), separated by a distance d = WOkpc. From equation 
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the nearby massive galaxy generates a gravitational acceleration on the order of 

\V(j)g\ ~ ^GMgao/d ~ 5.36 -1 ao, (16) 

within the — m system. The period of the test particle is appreciably less than that of the 
dwarf galaxy's orbit. Possible resonances between the dwarf galaxy's orbital period and that of the 
test particle are neglected. According to equation (|9|), equation (||) is valid for distances greater 
than r m i n = 1189pc and we choose 

r « 2378pc. (17) 

Tidal accelerations for distances such as r are neglected, so that V0 9 is considered to be a constant 
external gravitational acceleration within the — m system. The geometry is shown in Figure |l|. 

If the initial conditions are chosen as z = p z = 0, according to equation ([?]), Kepler's Theory 
is recovered with an effective gravitational constant (eq. @). These initial conditions correspond 
to orbits constrained to the equator of and are known as equatorial orbits. (The equatorial 
plane is defined as the plane containing M^, perpendicular to the z direction (see Fig. [I]).) It then 
follows that 

l 2 z = /x r Aa(l-e 2 ) (18) 
H _ IS—*, (19) 

T = 2 m 3 '\f^, (20) 
V A 

where a is the semi-major axis, e the eccentricity, T the period and E, the energy. In Newtonian 
theory, the relation between the period and the semi-major axis is Kepler's third law. The 
equatorial orbits are uniquely defined by equations (|l8|) and (|l9|). 

In the MOND limit, the period predicted by equation (||), for equatorial orbits, follows from 
equations ® and | 



27rr 3 / 2 i 



G(M d + m) : 



which, according to equations (15), (16) and (|17|), results in 

T ss 5.6 x 10 8 yr. (21) 



The physical configuration defined by equations (|T^) and ([0]) is described in the following. 
We assume the MOND limit in equations (||) and (||) and choose a g = 7/16, corresponding to the 
linearized lower order resonance of 3/4 in Table [p. Under these conditions, the surfaces of section 
for the phase orbits of equation (0) (Figs. || - 1|) are obtained for an energy E = —6.817 x 10 50 erg, 
which corresponds to distances of the order of that in equation (|l~7|). As a numerical check, the 
energy was found to be conserved to one part in 10 10 . 
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Regions of instability increase for smaller values of l z . Equatorial orbits are denned as 
z = p z = in equation (^). For equatorial orbits, l z is connected to the eccentricity by equation 
(|i~8p. According to equation (|14|), for equatorial orbits with sufficiently small eccentricities, almost 
the entire phase space is occupied by KAM tori surrounding the stable orbit, as shown in Fig. ^. 
The surface of section in Fig. [2] is obtained for l z f fi r = 57.09kpc km/s, which corresponds to an 
equatorial orbit with eccentricity e = 0.4. It can be seen in Fig. || that this particular equatorial 
orbit is stable. 

The MOND potential given in equation (|3|) is valid for distances greater than r m j n , defined in 
equation (||). This potential predicts equatorial orbits which belong to a plane perpendicular to 
the external gravitational acceleration produced by M g (see Fig. |I]). We numerically found that 
the equatorial orbits are stable for eccentricities e < 0.65. 

With decreasing values of l z , for example, l z /fi r = 41.20kpckm/s, the equatorial orbit becomes 
unstable. There exist other stable periodic orbits, which are not in the plane z = (see Fig. |3|). 

Associated with the upper and lower stable islands in Fig. ||, there is an exact 2 : 1 resonance, 
which is a periodic orbit in the meridional plane, defined by p, z coordinates. In three dimensional 
space x, y, z, this periodic orbit corresponds to a quasi-periodic orbit in a spatial region, which 
is a rounded cylindrical shell surrounding M^, as shown in Figure ||. A closed orbit in three 
dimensional space requires an exact resonance between the azimuthal, vertical and epicycle 
frequencies. 



4. Conclusions 

There exists a previously calculated potential ip in the MOND theory, given by equation @, 
for a free falling sphere of mass in a constant external gravitational acceleration V0 9 . p is 
valid when |V^» S | >> |V</?|. We assume that the existence of a second mass m « does not 
modify the potential. 

The system — m, bounded by the potential tp, is described by the Hamiltonian given 
in equation ([?]). Since the potential is axially symmetric, only the component of the angular 
momentum in the direction of V(/> s , l z , is conserved. As a consequence, when z = p z = and l z 
is the total angular momentum, motion is periodic and occurs in a plane perpendicular to V</> 5 , 
which we call the equatorial plane (Fig. |l]). Kepler's theory is then recovered with an effective 
gravitational constant (eq. (pi)), and l z is related to the eccentricity by equation (18). 



We analyze the phase space of H, using the numerical technique of surfaces of section. It 
is found that for eccentricities e < 0.65, the equatorial orbits are stable with respect to small 
variations in the initial conditions (Fig. ||). For smaller values of l z , the stability of the equatorial 
orbit is lost and there is an increase in the size of the chaotic regions (Fig. ||). There is a disk 
perpendicular to V</> s , as well as some regions not in the disk, which are stable against variations 
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of initial conditions. 
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Table 1. Resonances 



7/16 3/4 

9/25 4/5 

11/36 5/6 

24/49 5/7 



-10- 



REFERENCES 

Arnold, V.I. 1989, Mathematical Methods of Classical Mechanics, (New York: Springer- Verlag) 
Arnold, V.I. 1961, Soviet Math. Dokl., 2, 501 
Arnold, V.I. 1963, Russian Math. Surv., 18:6,85 
Bekenstein, J. D., Milgrom, M. 1984, ApJ, 286, 7 

Bekenstein, J. D. 1987, in: 2nd Canadian Conference on General Relativity and Relativistic 
Astrophysics, ed. C. Dyer, (Singapore: World Scientific) 

de Block, W.J., McGaugh, S.S. 1998, ApJ, 508, 132 

Brada, R., Milgrom, M. 1999, ApJ, 519, 590 

Danby J. M. A. 1971, Qualitative Methods in Celestial Mechanics, (Sao Paulo: Publicacoes do 
Instituto de Matematica e Estati'stica da Universidade de Sao Paulo) 

Faber, S.M., Jackson, R.E. 1976, ApJ, 204, 668 

Henon M., Heiles C. 1964, AJ, 69, 73 

Kolmogorov, A.N. 1954, Dokl. Akad. Nauk. SSSR, 98, 527 

Lichtenberg A.J. & Lieberman M.A. 1983, Regular and Stochastic Motion Second Edition, (New 
York: Springer- Verlag) 

Lindstedt, M. 1882, Astrom. Nach., 103, 211 

McGaugh, S.S., de Block, W.J.G. 1998, ApJ, 499, 66 

Milgrom, M. 1983a, ApJ, 270, 365 

Milgrom, M. 1983b, ApJ, 270, 371 

Milgrom, M. 1983c, ApJ, 270, 384 

Milgrom, M. 1989, ApJ, 338, 121 

Moser, J. 1962, On Invariants Curves of Area-Preserving Mappings on an Annulus, Nachr. Akad. 
Wiss. Gotingen, Math. Phys. Kl, p.l 

Moser, J., Siegel, C. 1971, Lectures in Celestial Mechanics, (Berlin: Springer- Verlag) 

Oort, J.H. 1963, The Moon, Meteorites, and Comets, ed. B. M. Middlehurst and G. P. Kuiper 
(Chicago: University of Chicago Press) 

Ostriker, J.P., Peebles, P.J.E. 1973, ApJ, 186, 467 



- 11 - 

Salcedo, F.J.S., Gamez, A.M.H. 1999, A&A, 345, 36 
Tully, R.B., Fisher, J.R. 1977, A&A, 54, 661 



This preprint was prepared with the AAS IAT^X macros v4.0. 



-12- 



Fig. 1. — Schematic diagram, showing the two masses = 1O 8 M0 (e.g., nucleus of a dwarf galaxy) 
and m = 1O 5 M (e.g., a globular cluster) and the neighboring galaxy M g = 5 x 10 11 Mq (e.g., the 
Milky Way). According to the potential given in equation (|2j), the M^ — vn system is gravitationally 
bounded to M g , at an approximate distance of d = lOOkpc. The external gravitational acceleration 
field due to M g is \V(f> g \ ~ 5.36 _1 a , where a = 2 x l(T 8 cm/s 2 is the MOND constant. The z 
direction is the projection of the relative distance r between and m in the direction of the 
external field V0 g and p is the orthogonal projection. The equatorial plane is defined by z = and 
is shown in gray. r m j n is the minimum distance for which the potential in equation (|3|) is valid. 

Fig. 2. — Surface of section, for E = —6.817 x 10 50 erg, and l z /p r = 57.090kpc km/s. 

Fig. 3. — Surface of section, for E = —6.817 x 10 50 erg, and l z /fA r = 41.201kpc km/s. 

Fig. 4. — Orbit in (x, y, z) space in units of kpc, corresponding to the 2 : 1 resonance of Figure ||. 



